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Abstract: The generation of THz radiation from ionizing two-color fem- 
tosecond pulses propagating in metallic hollow waveguides filled with Ar is 
numerically studied. We observe a strong reshaping of the low-frequency 
part of the spectrum. Namely, after several millimeters of propagation 
the spectrum is extended from hundreds of GHz up to ~ 150 THz. For 
longer propagation distances, nearly single-cycle near-infrared pulses with 
wavelengths around 4.5 jj.m are obtained by appropriate spectral filtering, 
with an efficiency of up to 0.25 %. 
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1. Introduction 

In recent years, the range of wavelengths where coherent radiation can be generated has grown 
dramatically into both high and low frequency domain. Remarkably, most of the methods to 
obtain radiation at extreme frequencies use, in one or the other way, nonlinear processes in 
laser-induced plasma. One prominent example is high harmonic generation (HHG) where fre- 
quencies thousand times larger than the frequency of the pump pulse are excited, exploiting 
the recollision dynamics of electrons ionized by the intense light pulses. More recently it was 
demonstrated that a two-color fs beam allows generation of new frequencies just in the oppo- 
site part of the spectrum, namely in the THz range, hundreds times smaller than the optical 
pump frequency. To this end, a short two-color pulse of fundamental frequency and second 
harmonic is strongly focused into a plasma spot |[T HT2ll The observed THz emission generated 
in this scheme has been attributed to the laser-induced plasma current in the asymmetric two- 
color field [7.8]. Using such scheme generation of strong THz radiation was reported, with a 
spectrum which can be as broad as 70 THz |8). Such broad-band coherent radiation can allow 
new applications providing the possibility to probe complex molecules or as an analytical and 
imaging tool in a broad range of fields. 

In this article we consider a modification of the above described focusing geometry in a bulk 
gas by using a metallic hollow waveguide with a cladding from aluminum to guide both THz 
and optical radiation. The use of a waveguide prevents the diffraction of light for all wave- 
lengths involved. In addition, such setup allows to realize nearly single-mode operation for a 
wide range of frequencies. High threshold intensities for ionization-induced processes in the 
range larger than 100 TW/cm 2 have been already realized in the context of HHG schemes in 
hollow waveguides II13II14I . We will show that during propagation a dramatic spectral broad- 
ening of the low-frequency part of the spectrum occurs, caused by only modest changes in the 
fundamental and second harmonic fields. In particular, ultrashort pulses of 10 fs duration lead 



to the formation of low-frequency spectrum extended from several hundreds of GHz to approx- 
imately 150 THz. Upon further propagation we report the generation of nearly single-cycle 
pulses with central wavelengths around 4.5 [Xm and efficiency of the order of 0.25 %. 

2. The model 

In the previous works, THz emission by two-color fields in gases has been interpreted by four- 
wave mixing rectification 0}[5]]. In contrast, later studies attributed THz emission to the laser- 
induced plasma-current in the asymmetric two-color field ll7 Ul2l described by the local quasi- 
classical model of the electron current in the laser field without taking into account propagation 
effects. As an alternative approach, particle-in-cell simulations were used l9l [l0l . Those meth- 
ods are, however, computationally very expensive, and propagation effects have been taken 
into account over distances of a few micrometers only. In such models the microscopic cur- 
rent is described as a sum over contributions from all electrons, born at discrete times t n . The 
"continuous" analog of this approach is given by ( ll7l[8l [Tll ') 

t 

J(t) = q J v(t,t )p(t )dt , (1) 

where q is the electron charge. The current density J(?) at the time t is obtained by integrating 
over all impacts from different ionization times fo: \(t,to) is the velocity of electrons at time t 
which were born at time to, and 8p(to) = p(to + 8t ) — p(?o) ~p(to)St is the density of electrons 
born at time to Q. This approach was utilized in the previous works of Kim and coauthors 

GHTQ- 

The additional integration over the variable to increases the numerical effort, but one can 
eliminate the microscopic velocity distribution v(?,?o) and derive a differential equation for 
the macroscopic current density J(f). Assuming no interaction with other particles and zero 
velocity of new-born electrons, the equation for the electron velocity can be written as v(f ,?o) = 
J- LE(t)c/t (where E is the electric field, m e is the electron mass). With the simple additive 
property of integrals we can rewrite it as v(f ,fo) = v(f, — °°) — v(fo, — °°). Substituting this into 
Eq. (HJ and differentiating with respect to t we obtain J(r) = q\(t)p(t). It is easy to see that 
by introducing new variables v' = vexp(— V e t), J' = Jexp(— v e t ) (where v e is the electron-ion 
collision rate) and proceeding analogously we can include collisional effects as well, and finally 
obtain 

J(t) + vJ(f) = — E(t)p(t). (2) 
m e 

Indeed, Eq. (f2]) coincides with the well known (and very commonly used) equation for the 
current in plasmas H15II16I . 

In the present article, the static model for the tunneling ionization ifTTl is used: 

p(t)=W ST [E(t)][p, t -p(t)}, w sr [£(r)] = a^yexp|-j3^J. (3) 
Here, p at is the neutral atomic density, and E a w 5.14 x 10 11 Vm~' the atomic field; the coeffi- 

5/2 3/2 

cients a = A-(O a r^ and [5 = (2/3)r H are defined through the ratio of ionization potentials of 
argon and hydrogen atoms r# = Uav/Uh, UAr = 15.6 eV, Uh = 13.6 eV, and co a « 4.13 x 10 16 
s- 1 003. 



We use the convention / = ^ . 



To describe the light propagation, we use the so called forward Maxwell equation (FME) 
for the fast oscillating optical field E(r,z,f) in the waveguide with the optical axis along z 
direction JT8): 
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where r = {x,>'}, A± = d xx + d yy , n((o) is the refractive index of Ar, k(co) = n{(0)(o/c. E(r,z, ft)) 
and P„/(r,z, ft)) are the Fourier transforms of E(r, z,f), P„/(r,z,f) with respect to time, P„/ is 
the nonlinear polarization, which includes both the Kerr (Pfcerr = £oX^ |E| 2 E) nonlinearity and 
the plasma contribution: P„; = pKerr + 'J/ to + iJi OS sl to, where J/ OSJ = Wsx(p a t — P)Ua,-/E is a 
loss term accounting for photon absorption during ionization; "division by a vector" means the 
component-wise division 1171 . Decomposing the field into a series of linear eigenmodes of the 
waveguide {F 7 (r), j = 1, ... °o} with corresponding propagation constants j3/(fi>), we obtain 
the following set of equations for the amplitudes of the eigenmodes Ej, coupled through the 

nonlinear polarization term P~!' = f FyP„/dr: 



3. Numerical simulations 

In our simulations, we assume a dielectric waveguide with aluminum coating of the inner walls 
and diameter d = 100 fim. Such waveguides can be routinely produced (see, e.g., lfl9l ). The 
main purpose of the waveguide is not to modify the dispersion relation but to confine both 
pump and THz radiation (especially the latter, because it quickly spreads out otherwise as a 
result of diffraction [25]). Taking into account that THz as well as fundamental and second 
harmonic frequencies undergo strong absorption due to plasma generation a long waveguide is 
not useful. Therefore, we restrict our simulations to distances < 1 cm. 

The linearly -polarized waveguide mode EHn has no cut-off at low frequencies, and the fun- 
damental and second-harmonic frequencies can be coupled into the waveguide with almost 
100 % efficiency from the free-space propagating laser mode. The dispersion and losses for 
the mode EHn of the 100-jUm waveguide are shown in Fig. [TJb) in low-frequency range. 
Both quantities are calculated using the direct solution of the corresponding boundary value 
problem ||20l . Drude model for the dispersion in aluminum ETl as well as Sellmeier formula 
for the dispersion in argon |22l . One can see that the losses are relatively large only in the 
small frequency range corresponding to wavelengths of the order of the waveguide diameter. 
The aluminum coating was chosen because then the modeling of linear dispersion and losses 
in the whole frequency range from sub-THz to sub-PHz is straight forward using the Drude 
model I2TI . In contrast, for dielectric coatings (such as fused silica) such simple description for 
the whole spectral region of interest fails. Nevertheless, we expect even weaker losses than for 
the metallic waveguide. 

Results from simulations of ultrashort pulses with 10 fs duration and frequencies Vo = 
375 THz (corresponding to the wavelength A = 800 nm) and 2vo , having intensities I a = 
10 14 W/cm 2 and ha = 2x 10 13 W/cm 2 correspondingly are shown in Fig. |2] We checked 
that the transverse field distribution consists predominantly of the fundamental and few higher 
order modes. At the initial stage of propagation, broadband THz radiation is generated from 
approximately 0.5 THz to 80 THz (see Fig. |2ta) and (b)). With further propagation (up to 
2 mm) the low-frequency part of the spectrum grows in the direction to higher frequencies (up 
to 150 THz). At the same time, the spectral intensity around 10 THz degrades because of the 
strong losses of the waveguide mode. From 2 mm up to ~ 7 mm the long-wavelength part 
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Fig. 1. (Color online), (a) Scheme for THz generation considered in this article. A short 
pulse with center frequency Vo = 375 THz (corresponding to wavelength A = 800 nm) 
and its second harmonic are focused into an aluminum coated hollow waveguide. The THz 
emission is detected at the exit of the waveguide, (b) Losses a (red solid line) and refrac- 
tive index change 8n = kc/co (green dashed line) induced by the Al waveguide with the 
diameter d = 100 jlm for the mode EHu versus frequency fl = (aj1%. 
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Fig. 2. (Color online). Evolution of the THz spectrum upon propagation for a waveguide 
with diameter d = 100 flm, input intensities I m = 10 14 W/cm 2 , h m = 2 x 10 13 W/cm 2 
and pulse duration 10 fs. (a) The spectrum and (b) low-frequency part of spectrum (note 
the logarithmic scale of y-axis in (a) and both axes in (b)) for z = 0.5, 2, and 10 mm 
of propagation (black dashed, green dotted and red solid lines), (c) Pulses obtained by 
filtering the long-wavelength part of the spectrum below 250 THz. Propagation distances, 
color coding and line styles are the same as in (a). 



of the spectrum shifts continuously towards higher frequencies. Finally, starting from 7 mm 
the spectrum does not change considerably anymore, except for its decrease due to losses. At 
this propagation distance, the long-wavelength part of the spectrum is localized around 60-80 
THz, with the center around 66 THz (corresponding to a wavelength ~ 4.5 jxm). The observed 
efficiency of this "frequency down-conversion" reaches 0.25 percent. 

The electric field obtained by filtering the long-wavelength part of the spectrum for different 
propagation distances is shown in Fig.|2c). One can see that despite the strong reshaping of the 
spectrum, the envelope of the obtained THz pulses does not change significantly. The amplitude 
of the pulse grows considerably during propagation, from approximately 1 % of the pump 
amplitude at z = 0.5 mm (which correspond to the conversion efficiency ~ 10~ 4 ) to about 5 % 
at z = 10 mm (which correspond to the conversion efficiency ~ 2.5 x 10~ 3 ). 

Fig. |2|a) reveals that in contrast to the strong reshaping of the low frequency components, 
the spectrum in the range of the pump frequencies remains almost unchanged. Nevertheless, the 
spectrum of the pump broadens slightly upon propagation. Due to plasma-induced refraction 
index change the maximum frequency of the fundamental pulse is slightly shifted to the blue 



123II24I . In our case the shift is small, only about 10 THz. As elaborated in 1251 . even small 
changes in the pump fields can affect dramatically the spectral shape of the generated low- 
frequency radiation. Such changes affect significantly the temporal positions of the ionization 
events. Contributions from the electrons born in different ionization events are superimposed, 
and interference effects lead to a reshaping the THz spectra. 




200 400 600 800 1000 50 100 150 200 

Q. (THz) a (THz) 



Fig. 3. (Color online). Spectra for different input pulse durations (25, 50 and 1000 fs are 
shown by red solid, green dashed and blue dotted lines, respectively) for 1 cm of propaga- 
tion. In (a) and (b) the full spectra and their low-frequency parts are shown, respectively. 
Spectra are normalized to their maximum value, so that efficiencies can be compared. 

The situation is qualitatively similar for other pulse durations (see Fig. [3). In general, the 
efficiency of the THz generation is larger for shorter pulse durations (in Fig. [3] spectra are 
normalized to unity at their maxima, so the efficiencies can be directly compared). This obser- 
vation can be explained by the additional asymmetry in the current which is introduced by the 
finite pulse duration and therefore more pronounced for short pulses. 

4. Conclusion 

In conclusion, we studied THz generation by plasma currents through noble gas ionization in 
strong asymmetric laser fields. We have shown that the mechanism allows the relatively simple 
generation of extremely wide range of frequencies, starting from THz to near infrared. By using 
hollow waveguides we overcome the strong diffraction of the generated low-frequency gener- 
ation. We report that depending on the length of the waveguide either THz supercontinuum or 
almost single-cycle near-infrared pulses can be produced at output with ~ 0.25 % efficiency. 
We believe that our findings may open new intriguing possibilities to control the spectral and 
temporal shape of low-frequency radiation generated by the ionization dynamics in intense 
optical fields. 



